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Stellar archaeology: the evolving spectrum of FG Sge 
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Abstract. Over an interval of 120 years, the extraordinary object FG Sge has been transformed from a hot post-AGB star to a 
very luminous cool supergiant. Theoretically, this evolution has been associated with the reignition of a helium-shell during the 
post-AGB stage. A series of studies of the chemical composition of the photosphere have suggested that this evolution has been 
associated with a dramatic increase of approximately 3 dex in the abundances of s-process elements between about 1960 and 
1995. The problem with this apparent change is that it occurred at a time when the surface convection zone, which is governed 
by the star's effective temperature, could not have developed sufficiently to dredge processed material from the stellar interior 
to the surface. We have reviewed the chemical evolution of FG Sge by means of modelling the time- varying spectrum under a 
range of assumptions. By comparing these models with published observational data, a self-consistent picture has emerged. In 
particular we find that surface hydrogen has been depleted during the interval in question. In contrast, the s-process abundances 
have generally maintained a steady enhancement of around 1 — 2 dex, although some modest changes may have occurred since 
1950. This implies that FG Sge has not just completed dredging up freshly-produced s-process isotopes. However, there remains 
a contradiction between the observed hydrogen-depletion, the age of the associated planetary nebula, and current evolutionary 
models for a pre-white dwarf suffering a late thermal pulse. 

Key words, stars: fundamental parameters - stars: abundances - stars: evolution - stars: individual (FGSge) 



1. Introduction 

During the last century, FG Sagittae evolved from a hot cen- 
tral star of a planetary nebulae to become an F-type supergiant 
within approximately 90 years. Only two comparable objects 
are known; V4334 Sgr (Sakurai's object) covered a similar evo- 
lution in approximately five years to become, in 1999, a C-type 
giant smothered by its own dust shell (Duerbeck et al. 2002), 
and V605 Aql which, after a brief foray into the giant domain 
during the 1920's, returned to the high-temperature regime as a 
typical [WC]-type planetary nebula central star (Clayton & de 
Marco 1997). The rarity and brevity of such events arising from 
a cataclysm in the final evolution of low-mass stars is such that 
they are of profound interest. 

For these reasons, every observation of these stars is im- 
portant, not so much for their contemporary significance, but 
for their contribution to retrospective studies of evolutionary 
episodes which, even for such astronomically rapid events, may 
last longer than most scientists' careers. In particular, the devel- 
opment of increasingly sophisticated analytical methods render 
the original data of primary importance in any historical inter- 
pretation. 

The generally accepted interpretation for the redward evo- 
lution of FG Sge is that of a late thermal pulse in a young 
white dwarf or planetary nebula central star (Paczynski 1971; 
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Langer, Kraft & Anderson 1974; Sackmann & Despain 1974; 
Schonberner 1979; Iben 1984; Blocker & Schonberner 1997; 
Gonzalez et al. 1998; Lawlor & MacDonald 2003) wherein nu- 
clear helium-burning at the surface of the electron-degenerate 
carbon-oxygen core is reignited. Detailed models for these late 
thermal pulses have been generated by several investigators (cf. 
Schonberner 1979; Iben 1984; Blocker 1995), but the main 
features are principally the same. The energy released by the 
helium-shell flash forces the outer layers to expand in about 50 
to a few 100 years, depending on their mass. Hydrogen-rich, 
helium-rich and carbon-rich material is mixed by convective 
overshoot from the helium-shell flash, but the surface of the star 
retains its pre-pulse composition. A major component of this 
mixed region will be a supply of newly-produced s-process el- 
ements generated by the mixing of ^'^C and protons. As the star 
expands, its surface cools allowing surface convection to oc- 
cur. As the temperature bottoms-out, the convection zone digs 
deeper, dredging processed material to the stellar surface. 

The evolution of FG Sge in terms of effective tempera- 
ture ( Toff) and luminosity (L) has been established from long- 
term photometric studies (van Genderen & Gautschy 1995). 
These data appear to support the late thermal pulse hypoth- 
esis (Blocker & Schonberner 1997; Herwig 2001a; Lawlor 
& MacDonald 2003). Further evidence is the apparent en- 
hancement of s-process elements in the photosphere of FG Sge 
(Langer et al. 1974) as T^s dropped below 8 500 K, almost ex- 
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actly the point where surface convection is expected to develop 
(cf. van Genderen & Gautschy 1995). The dilemma, in sim- 
ple terms, is that the theoretical models for late thermal pulse 
evolution do not predict the dredge-up of s-processed mate- 
rial to the surface until some considerable time after the star 
has reached its minimum effective temperature. Moreover, the 
dredge-up should also be associated with a significant change 
in the abundances of hydrogen, helium, carbon and oxygen. 
Iben (1984) remarks that "The heretical conclusion to which 
one appears to be forced is that FG Sge has not just completed 
dredging up freshly produced s-process isotopes. If this is so, 
then the fault must lie in the estimates of surface abundances." 

Consequently, either the theoretical predictions or the ob- 
servations concerning the surface abundance of late thermal 
pulse products is wrong, or FG Sge is not a late thermal pulse 
product. 

More recently, a variant of the late thermal pulse phe- 
nomenon has been identified (Herwig et al. 1999)^ The "late 
thermal pulse" (LTP) model introduced above does qualita- 
tively well in reproducing the hydrogen-deficiency and high 
carbon abundances seen in hydrogen-deficient planetary neb- 
ulae central stars, but fails to reproduce the large oxygen abun- 
dances found therein. Herwig et al. (1999) found, however, 
that a "very late thermal pulse" (VLTP), occurring after the 
H-burning shell has become inactive, was able to produce a 
surface abudance pattern in good general agreement with that 
observed in the hydrogen-deficient post-AGB stars, primarily 
because he considered convective overshoot. The two models 
are different in one major respect. In the first, the convection 
zone produced in the He-rich shell at He ignition cannot extend 
into the hydrogen-rich envelope (Iben 1984). In the second, 
this convection zone is able to penetrate into the H-rich surface 
layers and transport protons downwards (Herwig et al. 1999). 
What is not immediately apparent from the published figures, 
although it is stated in the text, is that the nuclear-driven con- 
vection zones are able to propagate through to the surface be- 
fore the star cools sufficiently to develop a deep opacity-driven 
convection layer at its surface. This means that in the VLTP 
model, the disappearance of surface hydrogen and the appear- 
ance of highly-processed material can occur quickly and while 
the effective temperature is still significantly high. 

Another major difference between the LTP and VLTP mod- 
els is that the latter occurs on a much shorter timescale than 
the former. A VLTP occurs when the envelope mass is close 
to ~ 10~^ M0. Since the nuclear energy produced by protons 
ingested into the intershell region is enormous (^ 10'^'' erg in 
less than a week), expansion following a VLTP occurs on the 
very small thermal timescale (~ 1 y) of the tiny hydrogen- 
rich envelope (Herwig et al. 1997). In contrast, a LTP occurs 
when the envelope mass is also ~ 10"'* M©, but ~ 10~^ MQof 
hydrogen-free intershell must also be lifted, there is no proton- 
ingestion, and the expansion energy comes from the He-shell 
flash alone. Consequently the expansion occurs on a much 
longer timescale (^ 100 y). 



' In fact, both classes of thermal pulse were previously encountered, 
but not explicitly labelled, by Schonbemer (1979). 
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Fig. 1. The evolution of FG Sge in the L- T^s diagram (van Genderen 
& Gautschy 1995, Table 3: high Tes). A smooth curve has been drawn 
through the data. 

The VLTP model is therefore important for understanding 
the observed abundance evolution of FGSge and V4334Sgr 
(Gonzalez et al. 1998, Asplund et al. 1997). For example, it has 
been stated by several authors that, because of their respective 
evolution timescales, FGSge and V4334Sgr are examples of 
LTP and VLTP evolution respectively {cf. Herwig 2001a). 

It is the intent of this investigation to determine whether 
observation can indeed be reconciled with theory, and whether 
FGSge has suffered either a "late" or a "very late" thermal 
pulse. The next sections give a more detailed discussion of 
specific observations, followed by a presentation of the model 
atmospheres and spectra used and our subsequent analysis. It 
considers in more detail ideas introduced by Schonbemer & 
Jeffery (2003). 

The literature on FGSge is extensive, there are currently 
over 394 citations in the SIMBAD database between 1960 and 
the present (2006 July)^. It would be impossible and inappro- 
priate to survey critically every one of these articles, but inad- 
vertent omissions of material that others might consider semi- 
nal will no doubt have been made. However we recognise the 
contribution of all of these publications, and hence our prefa- 
tory remarks on the importance of primary data. These de- 
rive from the fact that much early spectroscopic material is no 
longer accessible to us. Consequently, any conclusions which 
run contrary to the original authors' are legitimately disputable. 
On the other hand, if it can be shown that a self-consistent pic- 
ture of FG Sge could be realised by a reassessment of these 
data, this paper will have achieved its goal. 

2. The expansion of FG Sge 

The evolution track derived from optical photometry and pul- 
sation analysis (van Genderen & Gautschy 1995, Fig. 1) is 
adopted as a principal benchmark. These authors actually give 
two evolution tracks, with low and high temperature scales re- 
spectively. We have adopted the high scale for reasons which 
will become clear later, although van Genderen & Gautschy 

^ There is only one citation in SIMBAD from before 1960: 
Hoffmeister's (1944) Ust of 171 new variables. 



C.S. Jeffery & D. Schonberner: Spectral evolution of FGSge 



3 



■ Her68 

♦ Lon74 

O stovg 

A AckSa 

A Mon90 

□ Gon98 

X TGG95 

o 7GG95 

• Jef03 



>o p 



4. 10 
4. 05 
4. 00 
3. 95 
3 90 
^ 3.85 

O 

3.80 
3. 75 
3. 70 
3. 65 
3. 60 

1955 1960 1965 1970 1975 1960 1985 1990 1995 2000 
Year 

Fig. 2. The run of T^g for FGSge as determined spectroscopically 
and photometrically. Both low ( x ) and high (dots) photometric Tcft 
scales are from van Genderen & Gautschy (1995: vGG95). The spec- 
troscopic measurements are from Herbig & Boyarchuk (1968: Her68), 
Langer et al. (1974: Lan74), Stone (1979: Sto79), Acker et al. (1982: 
Ack82), Montesinos et al. (1990: Mon90) and Gonzalez et al. (1998: 
Gon98). New results from the present study are also shown (JefOB). 
This and other figures appear in colour in the online version of this 
paper; parenthetical references to colour refer to the latter. 



(1995) prefer the low scale. Underlying this is the basic as- 
sumption that over a long timescale and in the absence of 
major events affecting the internal structure of the star, the 
evolution track will be smooth (Fig. 1). Superimposed on the 
evolution are short-term pulsational changes (AV^ ^ 0.15, 
A{B — V)^ 0.04), which will perturb instantaneous measure- 
ments of Toff and L. Other authors (Montesinos et al. 1990, 
Gonzalez et al. 1998) argue that the adopted B — V calibration 
is unreliable because of line-blanketing in the optical by en- 
hanced s-process elements. This does not apply to observations 
before approximately 1970, when the first s-process episode 
was claimed (Langer et al. 1974). In this paper we investi- 
gate the influence of s-process element enhancements on the 
{B — V) — Tcff calibration. 

It is also useful to check broad-band photometric measure- 
ments of TcB against other methods. Spectrophotometry and 
spectroscopy have both been used to measure Toff for FG Sge 
over the last forty years. 

Spectrophotometry of FGSge covering the decade 1968- 
1978 was calibrated against a sequence of MK standard su- 
pergiants (Stone 1979). Only wavelengths > 5 556 A were 
used in an effort to avoid the effect of blanketing by s-process 
elements. The effective temperatures inferred were ~ 500 K 
higher than those derived spectroscopically (Langer et al. 1974, 
Fig. 2). The systematic difference could be due to non-LTE ef- 
fects in the atmosphere or, more likely in the present authors' 
view, to the influence of rare-earth metal lines on the red spec- 
trum. These would have different consequences for the two 
methods of temperature measurement. However the measured 
cooling rate (310Kyr^^) agreed with contemporary measure- 
ments (SmoUnski, Climenhaga & Kipper 1976). This rate ap- 
peared to reduce later to 170Kyr~^ (Cohen & Phillips 
1980). 



We note that the spectroscopic measurements of effective 
temperature Toff do not always agree well with measurements 
of the ionization temperature Tjon measured from the same data 
(Langer et al. 1974) but are in satisfactory agreement with inde- 
pendent contemporary analyses (Chalonge, Divan & Mirzoyan 
1977). 

Spectral type alone can be a useful measure of Tcs if the 
calibration (e.g. Johnson 1966) is reliable, and gives results 
(Acker, Jaschek & Gleizes 1982) that agree very well with the 
photometric calibration (Fig. 2). It is important to note that the 
lowest temperatures are given for the interval 1975-1981. The 
authors are not aware of any subsequent spectroscopic mea- 
surement of Tcs that gives a value below 6 000 K. 

The addition of ultraviolet and infrared photometry during 
the 1980's provided more data with which to measure Tcff us- 
ing the total flux method. However the intrinsic short-term vari- 
ability of FG Sge made simultaneous UV, visual and IR pho- 
tometry difficult to achieve. The problem is most apparent in 
the infrared, where a 1 15 d intrinsic variation (Arkhipova 1993) 
was shown to be due to a 1 300 K fluctuation in Tcs and a factor 
two change in luminosity (Taranova 1987). 

The only study to exploit the lUE data (Montesinos et al. 
1990) used a simple comparison with standard model atmo- 
spheres (Kurucz 1979) to provide estimates of T^s, conclud- 
ing that Tcff = 6 500 K and log g = 1.5 in the interval 1982- 
1989. These authors also conclude that the previous expansion 
had stopped sometime in the early 1980's and that a possible 
rise in temperature had occurred. The questions raised by this 
study concern the appropriateness of the model atmospheres 
used, the influence of variability on the comparison of non- 
contemporary data and - as was to become clearer in the 1990's 
- the development of a dust shell masking the photosphere. 

Subsequently, colours in 1991-1993 were reported to be 
consistent with spectral type Gp, while the line spectrum 
remained essentially the same as in 1983 (Stone, Kraft & 
Prossner 1993). 1992 August marked the beginning of a new 
phase in the evolution of FG Sge (Papousek 1992) in which it 
has shown a series of deep obscurational minima reminiscent 
of an RCrB-type light curve (Royer 1999). The clouds of dust 
formed during these events have substantially transformed ev- 
ery subsequent observation of the underlying star, so that it has 
been difficult to follow any further evolution in, for example, 
Tcff or intrinsic brightness. 

Several authors have employed the intrinsic variability or 
pulsations of FGSge as a proxy for Tcs- The principle is 
that linear pulsation models provide a relation between mass, 
period, luminosity and temperature. Assuming that the star 
evolves at constant luminosity (and mass), then a pulsation pe- 
riod gives Tcff. However, the difficulty of establishing reliable 
periods and the necessary assumptions regarding mass and lu- 
minosity suggest these estimates should be used with caution. 

3. The chemical evolution of FG Sge 

The expansion of FG Sge has naturally caused huge changes in 
its spectrum. It has been argued that these are not purely due 
to changes in Tcff alone, but that there have been concomitant 
changes in the chemical composition of the photosphere. The 



4 



C.S. Jeffery & D. Schonbemer: Spectral evolution of FGSge 



most significant studies are reviewed here. It is appropriate to 
note that each study considers a different epoch, using a differ- 
ent observational dataset. Of equal note, each study approaches 
the abundance analysis using a different method with different 
underlying assumptions. The dangers of simply concatenating 
such a series of results must be self-evident. Attention is drawn 
to some of the major inconsistencies. 



3. 1. Herbig & Boyarchuk 1968 

The first data to permit a view of the chemical evolution 
of FGSge was provided by Coude spectrograms obtained 
at the Lick Observatory between 1960 and 1972 (Herbig & 
Boyarchuk 1968; Langer et al. 1974), together with some ear- 
lier low-dispersion plates from 1955. The dramatic spectral 
evolution from B4 1 to A5 la was initially interpreted as an ex- 
panding surface, or false photosphere, above an invisible cen- 
tral star (Herbig & Boyarchuk 1968). However, it is the ex- 
tensive tables of equivalent widths published by these authors 
that provide the most valuable and accessible data from which 
to reassess the history of FGSge. These include data from ab- 
sorption lines between 3300 and 4930 A, including those due 
to hydrogen and neutral heUum, Ught elements and iron-group 
ions. 

These equivalent widths were used by Herbig & Boyarchuk 
(1968) in a curve-of-growth analysis to measure abundances of 
20 species from hydrogen up to europium. The results were 
tabulated relative to a Cyg, after normalising all abundances so 
that the mean combined abundance of Si, Ti, Cr and Fe was 
the same for both FGSge and a Cyg (see appendix A). The 
following points are interesting: 

i) The absence of He i lines from 1962 onwards places hmits 
on Tes and helium abundance; 

ii) a probable excess of carbon (-1-0.7 dex) in most spectra; 

iii) a 0.5 dex overabundance of oxygen; 

iv) an apparent excess of cobalt and strontium (-1-0.8 dex) at all 
epochs; 

v) nearly aU other elements are overabundant by 0.2 to 0.3 dex. 

Assuming that the differential abundance measurements are 
relatively secure, a modern abundance analysis of a Cyg can 
be used (e.g. Albayrak 2000) to compute absolute abundances 
(see Appendix A). After correcting for substantial differences 
between a Cyg and the Sun, we could further conclude: 

i) carbon is overabundant by 1 dex, 

ii) helium, nitrogen - iron and nickel are near-normal, and 

iii) cobalt, strontium, zirconium, barium and europium are 
overabundant by from -hO.7 (cobalt) to -h2.1 dex (strontium). 

It is significant that large overabundances in heavy ele- 
ments, particularly cobalt and strontium, are evident from 1960 
onwards. We should, however, be cautious in evaluating these 
data: measurement errors were not given, and the curve-of- 
growth method has limitations. However, the equivalent widths 
themselves are invaluable as they permit further analysis using 
modem model-atmosphere techniques. 



3.2. Langer, Kraft & Anderson 1974 

As the star continued to evolve, it began to show abnormally 
strong Unes of Y II, Zr II, Ce II, La II and other s-process species 
sometime after 1967 (Langer et al. 1974). Apparent abun- 
dances rose to some 25 times (h-1.4 dex) the solar value. 
Commencing with a reanalysis of the final spectrum used by 
Herbig & Boyarchuk (1968), the most interesting results are: 

i) a drop to near-normal carbon abundance (low weight). 

ii) steady Ti, Cr and Fe abundances. 

iii) significant increases in Y and Zr. 

iv) the appearance of La, Ce, Pr, Nd and Sm with h-0.9 ... h- 1.5 
dex overabundances. 

It is difficult, however, to reproduce these values as the cru- 
cial equivalent widths have not been published, and only a pho- 
tographic reproduction of the coude spectra in a small region of 
the spectrum is accessible to us^. Taken together, these two pa- 
pers (Herbig & Boyarchuk 1968; Langer et al. 1974) pose a 
number of significant questions: 

i) What has happened to the carbon abundance in FG Sge? It 
appeared to start high and then fall. However, if dredge-up 
of s-process elements occurred during this epoch, then carbon 
should also have been enriched. This may not be a severe obser- 
vational constraint; the early carbon abundances were "based 
on one or two very high excitation Cl lines in an uncertain at- 
mosphere with kT ~ 0.5 eV, so the Boltzmann factor would 
be pretty sizeable, and the results pretty uncertain" (Kraft, pri- 
vate communication). Since the carbon abundance derived by 
Langer et al. (1974) is also quite uncertain, the case of the car- 
bon abundance as based on these two investigations remains 
unsettled. 

ii) Is FG Sge hydrogen-deficient and when did it become so? 
Is there any corollary with, say, V4334Sgr, claimed to show 
a similar transition to hydrogen-deficiency during its redward 
evolution (Asplund et al. 1997)? Deep envelope mixing and the 
enrichment of s-process material should be dominated by the 
dredge-up of helium and dilution of surface hydrogen. This is 
a very important question. A low hydrogen abundance would 
be difficult to detect after 1960, yet it would have a dramatic 
effect on the structure of the atmosphere and the strengths of 
many Unes in the spectrum. 

iii) Has FG Sge been s-process rich all the time? Are early en- 
hancements in Co and Sr symptomatic of a general s-process 
enhancement that was simply not recognised until the photo- 
sphere cooled sufficiently for ions with lower ionization poten- 
tials to become visible? 

iv) How would the adopted model atmosphere affect, for ex- 
ample, measurements of the change of abundance with effec- 
tive temperature? In the domain of interest, a very small drop 
in TeS produces a very sharp increase in W\ for the second 
ions of rare-earth elements. Would additional opacity from car- 



^ The authors have contacted two coauthors of Langer et al. (Kraft 
and Anderson) and the principal author of Cohen & Phillips (1980), 
who also had access to these data. Detailed records of the original 
work can no longer be located. Although the original plates should still 
be in the Lick plate archive, their extraction and analysis lie beyond 
our current resources. 



C.S. Jeffery & D. Schonbemer: Spectral evolution of FGSge 



5 



bon, or reduced opacity from hydrogen depletion be sufficient 
to cause a significant change in the original conclusions? 

3.3. Cohen & Phillips 1980 

Cohen & Phillips (1980) continued the work of Langer et al. 
(1974) using 4-m echelle spectra obtained each year from 1975 
to 1978. They introduced a ratio technique to deduce Tgff 
and chemical abundances from each spectrum. No lines with 
Wx > 200 mA changed by more than 40% during this epoch. 
They also found for strong lines that oc {Nij)'^, with 

/3 = 0.30 rather than the anticipated 0.50. Whether this is a 
consequence of the low surface gravity or of the composition 
of the photosphere remains to be verified. The authors' infer- 
ence from the upper limit to changes in W\ was that, in the 
absence of any coohng, the abundances of all elements with 
strong lines, including the rare earths, did not increase by more 
than 0.50 dex over the 3 yr time span. Any cooling would re- 
duce this maximum abundance change. One conmient is no- 
table; "These models . . . predict for solar abundances of the 
non-s-process elements (believed to be valid for at least the 
even iron-peak elements in FGSge) values of W\ smaller than 
those observed in the spectra." The authors suggest that this 
is due to the extended nature of the atmosphere. It is also the 
classic symptom of an atmosphere with a lower than usual con- 
tinuum opacity as would be occasioned, for instance, by a low 
hydrogen abundance. If it is not an abundance effect, the cause 
should have the same consequence for all fines, including the s- 
process elements, accounting for some, if not all, of their appar- 
ent overabundance. Unfortunately, due to severe line crowding 
in the blue, the spectra used in this study were obtained long- 
ward of 5000A so that the only Balmer line observed was Ha. 

3.4. Acker, Jaschek & Gleizes 1982, Cowley, 
Jaschek & Acker 1985 

Following the work at Lick Observatory, the spectrum of 
FGSge during 1979 - 1980 suggested a remarkable new phe- 
nomenon. In addition to the increase in s-process elements 
(barium and the lanthanides), the abundances of iron-group el- 
ements were significantly reduced (Acker et al. 1982; Cowley 
et al. 1985). 

3.5. Kipper & Kipper 1993, Kipper & Klochkova 
2001 

As well as abundance estimates for magnesium, four iron-peak 
elements and 13 rare-earth elements. Kipper & Kipper (1993) 
report the detection of C2 in 1992. The deduced abundance of 
carbon is dependent entirely on the assumed Jeff (and presum- 
ably on the assumed hydrogen abundance) but appears to be 
enhanced, lying in the range -2.70(5 500 K) < logA^'(C) < 
— 1.6(6 500 K), where A''(C) represents the fractional abun- 
dance (by number of atoms) relative to the total of all ele- 
ments. Consequently, carbon is considered to have increased 
along with or after the s-process elements. These authors used 
log 5 = 1.0 in their analysis. 



High-resolution spectra taken in 2000 and analysed with 
the same model atmosphere (5 500 K/ log 51 = 1.0) gave es- 
sentially the same results (Kipper & Klochkova 2001). The 
oxygen abundance, based on one single [O i] line, resulted in 
N{C)/N{0) = 2.5. Note that the model parameters imply 
log L/M ~ 3.5 (solar units). 

Although a hydrogen-deficiency was not reported in either 
study, it should be noted that the abundances deduced are Ukely 
to be substantially smaller if hydrogen is significantly depleted. 

3.6. Gonzalez et al. 1998 

Gonzalez et al. (1998) obtained modern high-resolution high 
S/N spectra between 1992 and 1996. Their fine analysis consid- 
ered primarily a spectrum obtained in 1994. From Fe i and Fe 11 
lines they deduced T^g = 6 500 ± 400 K, log g = 2.0 ± 0.5, 
(Fig. 3) and a 3 dex enhancement of s-process elements. About 
1 dex of this enhancement must have occurred between 1992 
and 1994 (Kipper & Kipper 1993). 

As with the photometric analysis of Montesinos et al. 
(1990), the difficulty posed by this analysis is that it requires 
the star to have both heated - possibly as it returns towards 
higher temperatures - and - to yield the measured surface grav- 
ity - shrunk in radius by a factor of ~ 10 since the observations 
reported by Acker et al. (1982). Such a change would have pro- 
duced a ~ 1.7 dex reduction in luminosity, but no concomitant 
change in V (4 magnitudes) was observed over this interval. 

It is therefore puzzling that emphasis should have been 
given to the 3 dex abundance change, but not to the radius 
change. The effective temperature and abundance measure- 
ments are not unaffected by this contradiction, since both de- 
pend on the surface gravity and model atmosphere adopted. 
Reviewing our own work on atmospheric fine analyses we 
find that, in general, lower surface gravities give lower effec- 
tive temperatures for the same ionization or excitation equifib- 
rium with d log Teg /d log 5 ~ 0.06. Reducing the Gonzalez 
et al. (1998) surface gravity by the requisite amount to give 
log L/M « 4.0 yields T^s ~ 5 500K. 

Secondly, if the abundances are as far from normal as these 
authors suggest, the implications of reduced hydrogen opacity 
and increased metal opacity will be extreme, as Gonzalez et 
al. (1998) were aware. The structure of the adopted ATLAS9 
model atmospheres with normal hydrogen and metal abun- 
dances would have been quite inappropriate for modelling the 
spectrum of such a peculiar star. With few reliable temperature 
indicators available to Gonzalez et al. (1998), one wonders if 
the starting premise that the Montesinos et al. (1990) results 
were "probably the most accurate published temperature esti- 
mates since the early 1970's" was unfortunate. 

However, even if both effective temperature and surface 
gravity are wrong, the abundance measurements may not be so 
badly in error. If they are not, then FG Sge will have taught us 
something important about the late stages of stellar evolution. 
It is therefore both legitimate and necessary to make a thor- 
ough reappraisal of the chemical evolution of FG Sge, includ- 
ing those measurements made with the most modem equipment 
and methods. 
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Fig. 3. The run of Tcb and log g for FG Sge (van Genderen & 
Gautschy 1995; solid lines) assuming M — 0.8 M© (below) and 
0.6 Mq (above). Solid squares represent a selection of converged 
STERNE models, including the low-gravity limits. Several models with 
intermediate gravities were also computed. The straight line (online: 
red) indicates a sequence of models with log L/A'I = 4.0, for which 
theoretical colour indices, emergent spectra and equivalent widths 
were computed. The dashed line (online: blue) shows part of the 
redward loop from the evolutionary track of a 0.625 M© post-AGB 
model suffering from a late thermal pulse (see section 6.1). The circle 
(online: blue) shows the model adopted by Gonzalez et al. (1998) for 
their 1994 October 21 spectrum (Section 3.6). 



4. Model atmospheres and synthetic spectra 

In order to interpret the evolving spectrum of FG Sge it is nec- 
essary to simulate the emergent radiation observed through its 
history under a range of assumptions. The principal constraint 
on these is the evolutionary track derived from optical colours 
and pulsation periods (van Genderen & Gautschy 1995). This 
track is shown in Fig. 3, assuming masses of 0.6 and 0.8 Mq, 
respectively. 

It will be seen that these tracks lie close to the classical 
Eddington limit and that obtaining suitable model atmospheres 
is not trivial. In the current era, the strong radiation field and 
extended atmospheres generally encountered in A-type super- 
giants would be considered to require a treatment of the failure 
of local thermodynamic equilibrium (LTE) and the curvature 
of the stellar photosphere, in addition to a luminosity-related 
stellar wind. However, in the present case, the quality of the 
available data are insufficient to warrant such effort so that the 
adoption of line-blanketed plane parallel model atmospheres in 
LTE may be considered a satisfactory approximation. 

The program STERNE (Schonberner & Wolf 1974; Jeffery, 
Woolf & Pollacco 2001) can compute low-gravity model atmo- 
spheres over a wide range of Tcs and chemical composition. 
Its use provides control over the assumed physics and allows 
us to monitor the convergence of each model atmosphere. We 
note that the current opacity distribution functions are incom- 
plete, being based on the Kurucz & Peytremann (1975) line list, 
and that the contribution of molecules to the equation of state 
and opacity in low T^s atmospheres is not considered. Since 
the calculations we make are largely comparative, the fact that 



they are used consistently to discuss all of the data for FG Sge 
across a large range in T^s far outweighs these limitations. 

Model atmospheres were computed assuming a standard 
(solar) mix of elements heavier than helium. Three grids were 
computed assuming number abundances of hydrogen and he- 
lium (nH,nHc) = (0.9,0.10), (0.10,0.90) and (0.01,0.99), 
which we have labelled hydrogen-normal, hydrogen-poor and 
hydrogen-deficient, respectively. The goal was to obtain a se- 
quence of models with log(L/ Lo)/(Af/ M©) = logL/Af = 
4.0, close to the observed value and spanning the range in Tcs 
exhibited by FG Sge since 1960. The grid computed (although 
not every model) is illustrated in terms of Toff and log g in 
Fig. 3. We note that FGSge's luminosity-to-mass ratio may 
even be larger than the one assumed here. 

Within each of these model grids, the composition 
of the heavy element component assumed in the calcula- 
tion of the equation of state and the continuous (bound- 
free) opacity was (by number fraction relative to the to- 
tal) : C:0.000363, N:0.000112, 0:0.000832, Mg:0.0000339, 
Al:0.00000269, Si:0.0000316, S:0.0000158, Ca:0.00000229, 
and Fe:0. 0000288, these being the standard solar abundances 
of these elements at the time they were adopted into sterne. 
The abundance of the most abundant element (H or He) was 
reduced by the total of these to ensure a total number den- 
sity of unity. A complete review of the opacity sources used in 
STERNE has been given recently by Behara & Jeffery (2006). 
The abundances used in the calculation of the line-opacity 
were those hard-wired into the opacity distribution functions 
(ODFs) - being the Kurucz (1970) ATLAS6 "pOO" ODE for the 
hydrogen-rich grids and the MoUer (1990) "he90" ODE for the 
helium-rich grids. 

Attempts to compute corresponding models using Kurucz' 
program ATLAS9 (Kurucz 1991) were not so successful. These 
generally failed at very low gravities due to negative pres- 
sures in the equation of state, although we note that Lamers 
& Fitzpatrick (1988) had computed much lower gravity mod- 
els using an older version of ATLAS (Kurucz 1979), achieving 
log .9 = 1.2 rather than 2.0 at T^ft = 10 000 K. 

Using model atmospheres from STERNE, the formal solu- 
tion program SPECTRUM (Jeffery et al. 2001) permits the com- 
putation of synthetic spectra (in LTE) over large wavelength in- 
tervals, the computation of individual line profiles and curves 
of growth and the computation of elemental line abundances 
for given equivalent widths. 

A database for transitions of light elements in hot stars is 
maintained by one of us (LTE_LINES, Jeffery et al. 2001). For 
the present study, substantially more extensive data for iron- 
group and heavier atoms are required. The iron-group data have 
been collated from the extensive line lists of Kurucz (1994), 
providing a total of over 144500 metallic absorption lines with 
wavelengths between 2900 and 8000 A. For s-process ele- 
ments, data for over 1 1 500 lines have been recovered from the 
Vienna Atomic Line Database (Piskunov et al. 1995; Kupka et 
al. 1999). Since this represents only a subset of the s-process 
element lines observed in FG Sge, the effects of these elements 
will be underestimated. 
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Fig. 4. Simulated colour indices (B ~ V) for a sequence of model at- 
mospheres with log L/M — 4 and normal hydrogen (solid: black), 
10% hydrogen (by number, solid: red), 1% hydrogen (solid: green) 
and normal hydrogen with enhanced s-process elements (broken: 
black). 

5. Key questions 

This section starts to address some of the key questions raised 
in the previous discussion. For example: 

- Why are there such substantial discrepancies in the T^s 
scales? 

- What was the evolution of the hydrogen abundance? 

- Is there any evidence of carbon enrichment ? 

- Did s-process enhancement ever really occur? 

5. 1. Are colour indices sensitive to composition? 

The STERNE model atmospheres include a calculation of the 
emergent spectrum sampled at 342 wavelengths from the X- 
ray to the infrared, with sampling approximately every 20 A 
in the UV and visible. A consequence immediately apparent 
is that for T^s < 10 000 K, the hydrogen-deficient models are 
substantially redder, in the optical region, than the hydrogen- 
normal models. 

These flux distributions have been augmented by com- 
puting detailed synthetic spectra with SPECTRUM using the 
linelists described above and 200,000 frequency points be- 
tween 2900 and 8000 A. These were convolved with appropri- 
ate detector response and filter transmission functions to com- 
pute theoretical intrinsic colour indices {U — B)o and {B — V)o. 
for models with both normal and modified abundances. In ad- 
dition to the variations in hydrogen abundances, we also com- 
puted colours for a models in which the abundances of s- 
process elements strontium, yttrium, zirconium, barium, lan- 
thanum, cerium, praseodymium, neodymium. samarium, eu- 
ropium, terbium, dysprosium, erbium, ytterbium, lutetium, 
hafnium, and lead were enhanced by two dex. Most of these 
have been reported at some time to be strongly overabundant 
in FG Sge. 

The results for {B — V) are shown in Fig. 4 which confirms 
that the cooler hydrogen-deficient models may be redder at a 
given temperature than their hydrogen-rich counterparts. This 
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Fig. 5. Simulated colour indices (17 — B) vs. {B — V) for a sequence 
of model atmospheres with log L/M = 4.0: 1) normal hydrogen 
(black), 2) 10% hydrogen (by number, red), 3) 1% hydrogen (blue) 
and 4) normal hydrogen with s-process elements increased by 2 dex 
(black dotted). The hydrogen number fractions are marked at the cool 
end of each sequence. The effective temperatures of selected models 
are labelled on the 1% hydrogen sequence, models with the same Tes 
are marked in each sequence. The evolution of FGSge from 1960 to 
1985 is also shown (van Genderen & Gautschy 1995, broken black). 
The standard galactic reddening is shown as an arrow. 



is a direct consequence of the much lower continuum opaci- 
ties in the hydrogen-poor models. Strong metal line blanketing 
due to singly-ionized iron-group elements in the near-UV co- 
incides with the peak of the Planck function and has a much 
stronger effect than in the hydrogen-rich model atmospheres - 
the ratio of line to continuous opacity is the crucial quantity. 
The increase in metal-line blanketing is partially compounded 
in the hotter models by a reduction in the contribution from the 
high-order Balmer series and Balmer continuum. One impli- 
cation is that, for observations of FG Sge, the assumption of a 
normal hydrogen abundance could lead to an underestimate of 
Teff by as much as 1 000 K for T^s < 9 000 K. 

The original object of this experiment was to check the ef- 
fect of enhanced s-process element abundances on the colour- 
Tcff calibration. Although only a subset of all s-process lines 
was included. Fig. 4 shows that these effects are probably not 
severe, contrary to the assertion by Montesinos et al. 1990. 
They produce much smaller changes in (i? — V^) than caused by 
changes in the hydrogen abundance. Therefore, if the hydrogen 
abundance is normal, the existing Tes — {B — V) calibration 
(van Genderen & Gautschy 1995) should be acceptable. 

5.2. Can pfiotometry solve the hydrogen problem? 

Taking the previous discussion one step further, it is interesting 
to ask whether photometry alone can constrain the chemical 
evolution of FG Sge. It has been shown that the Teff scale de- 
pends on the actual hydrogen abundance. For early-type stars, 
the Balmer continuum dominates the position of a star in the 
two-colour iU — B) — [B — V) diagram. Figure 5 compares 
the locus of hydrogen-rich, poor and deficient model sequences 
with the evolutionary track of FG Sge. 
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While reducing the hydrogen from near-normal to ~ 10% 
reddens the evolutionary track in (B — V), as discussed above, 
the concomitant reduction in the Balmer continuum makes the 
tracks in the two-colour diagram distinguishable only in se- 
lected temperature and abundance domains. 

The only clear conclusion is that, assuming a standard ex- 
tinction law, FG Sge has been too red for a normal composition 
since at least 1975. However, since the star has been losing 
mass, its continued reddening may have as much to do with 
nonstandard extinction as with changing composition. 

5.3. How robust is the photometric calibration? 

If Tgff is inferred from [B — V) alone, one has to compute 
{B — V)o assuming some value for the reddening E^b-v- 
Most authors have adopted Eb-v = 0.4, apparently following 
Herbig & Boyarchuk (1968) who compared observed colours 
with those expected for the current spectral type. While these 
authors also reported an unusual value for E'u-b/ Eb-v = 
0.4 rather than the more usual value between 0.7 and 0.8, there 
have been no other suggestions of anomalous extinction before 
1990, when dust-ejection episodes became frequent. It will be 
seen that adopting .Eb-v = 0.4 is quite consistent with the op- 
tical/UV flux distribution in the 1980's and discussed in Section 
5.7. 

5.4. Was the composition of FGSge unusual before 
1965? 

The first step toward determining whether a substantial change 
in surface composition actually occurred is to establish a zero 
point. Specifically, was already FGSge depleted in hydrogen 
or rich in carbon or s-process elements in 1960-1965? 

Since Herbig & Boyarchuk (1968) provide equivalent 
width measurements for over 300 lines and five different spec- 
trograms, it is possible to assess this question quantitatively. 
Rather than attempt a detailed reanalysis of each spectrum, it 
is preferable to establish what systematic effects can be estab- 
lished. 

We have computed predicted equivalent widths for 278 of 
the lines measured by Herbig & Boyarchuk (1968) for model 
atmospheres with log L/M = 4.0 in the range 14 000 > 
Tcff > 5 000 K and with hydrogen compositions of 90%, 10% 
and 1% by number. The abundances of all species other than 
hydrogen and helium were assumed to be solar. We also com- 
puted a set of equivalent widths for the 90% hydrogen models 
in which all species other than hydrogen and helium were en- 
hanced by a factor of ten. Microturbulent velocities Vt = 5 and 
10km s~^ were used. In Figs. 6 and 7, the equivalent widths 
for selected fines are plotted as a function of temperature, with 
Vt = 5kms~^. The observed equivalent widths for each fine 
are also plotted, assuming T^s given by the photometric cal- 
ibration (van Genderen & Gautschy 1995, Fig. 1). The maxi- 
mum uncertainty in T^ff is ±500 K. In particular, these figures 
demonstrate how sensitive different lines are to changes in Tgs, 
hydrogen abundance and the elemental abundance. 



Table 1. Elemental abundances used in Fig. 8, written as log N^i — 
logiVH -t- 12, compared with solar photospheric abundances by 
Asplund, Grevesse & Sauval (2005). 



Elem 


FGSge 


Sun 


Fe 


7.25 


7.45 


Sr 


3.78 


2.92 


Y 


3.56 


2.21 


Zr 


3.88 


2.59 


Ba 


2.80 


2.17 


La 


2.69 


1.13 


Ce 


2.85 


1.58 


Sm 


2.00 


1.01 



A recurring question through the following analyses is that 
of the microturbulent velocity. The appropriate value to choose 
for this artificial factor is related to T^g and log g (Gray, 
Graham & Hoyt 2001), and also to the specific ion in question 
(Albayrak 2000). To constrain Vtina problem where the num- 
ber of degrees of freedom is comparable with the quantity of 
available data, we have chosen vt to be either 5 or 10 kms~^ 
according to which gives the most conservative interpretation. 
In practice, this meant using vt = 5 km s~ ^ when FG Sge was 
an early A star, and vt = lOkms"^ subsequently. 

A number of strong conclusions concerning the foUowing 
elements in FG Sge in 1960-1965 may be drawn: 

- He is normal. 

- C and O are not detectably different from normal, but C/O 
> 1 cannot be excluded. 

- Sr is significantly above normal (~ x 10). 

- Y and Zr are not significantly different from normal. 

Weaker conclusions include the following: 

- Mg, Si, and Ca are not well determined at these tempera- 
tures. 

- Sc and Ti appear normal. 

- Mn, Fe, Co and Ni may be somewhat underabundant. 

- Ba and Eu are significantly above normal (~ x 10, one ab- 
sorption fine only for each species. 

These conclusions are in line with the differential measure- 
ments given by Herbig & Boyarchuk (1968) if the abundance 
differences between a Cyg and the Sun are also taken into ac- 
count. 

5.5. Was there an s-process enhancement episode 
around 1965-1970? 

If some s-process elements were already enhanced in the in- 
terval 1960-1965, what is the significance of the dramatic in- 
crease reported between 1965 and 1970 (Langer et al. 1974)? 
The evidence rests on a reanalysis of the plate obtained in 1965 
by Herbig & Boyarchuk (1968) in which both fight and heavy 
s-process elements are reported to be low. It is appropriate to 
note the remark by Langer et al. that "In a first reconnaissance 
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Fig. 6. Equivalent widths of selected lines as a function of T^a for model atmospheres with log L/M = 4.0 and nn = 0.9 (solid), 0.1 
(dashed), and 0.01 (dot-dashed), assuming solar abundances and riHe = 1 — riH, and with nn = 0.9 and 10 x abundances (dotted, except 
riHe). Measured equivalent widths for FGSge during 1960-1965 are shown with symbols (Herbig & Boyarchuk 1968). 



of the abundance changes ... , rather crude methods seem en- 
tirely appropriate ... ". The 1965 spectrum only contains one 
line of barium amongst the heavy elements considered to be 
normal in 1965 by Langer et al. 

We have already found this line and one due to europium 
to be stronger than expected. Whilst we find that yttrium and 
zirconium are not detectably different from normal in 1965, we 
do find another light s-process element, strontium, to be about 
1 dex overabundant in 1965, before the apparent enhancement 
episode began. 

We have considered the question from two angles. First, we 
have calculated theoretical equivalent widths of several lines 
for the given s-process abundances and T^g. By deducing what 
Tgfi or Vt would be required if the abundances were to be nor- 
mal and yield the same equivalent widths, we conclude that it 
would have been be impossible to confuse an s-process over- 
abundance with a systematic error in T^s or a significant dif- 
ference in Vt between FG Sge and the comparison stars a Per 
and a Cyg. 

Second, we have attempted to reconstruct the photographic 
spectra presented by Langer et al. (1974) in Fig. 8. The loga- 
rithmic greyscale images represent the region 4398 to 4640 A, 
degraded to an instrumental resolution of 0.25 A judged to cor- 



respond approximately to 8 A mm~^ dispersion plates. Spectra 
were computed for Toff from 8 000 K to 5 500 K in steps of 
500 K, for two chemical mixtures, one with solar abundances 
and one with abundances reported for 1972 August (Langer 
et al. 1974, Table 1), and for vt = 5 and lOkms"^ Those 
with Vt = lOkms"^ are shown in Fig. 8. The temperature se- 
quence includes spectral classes from A8I to GO I, bracketing 
the Langer et al. range FO 1 to F6 1. The progressive enhance- 
ment with falling T^s of lines due to s-process elements com- 
pared with lines of iron-peak elements is evident, and confirms 
that s-process elements were overabundant. 

From the limited data available, we draw the following con- 
clusions for the interval 1965-1972. 

- the Langer et al. results are broadly reliable. 

- abundances of Ught and iron-group elements are near nor- 
mal. 

- Y and Zr are ~ 1 dex overabundant and had apparently 
increased after 1965. Note that Sr had been ~ 1 dex over- 
abundant before 1965, but was not measured subsequently. 

- Ba, La, Ce, Pr, Nd and Sm are ~ 1 dex overabundant; only 
Ba had been measured previously when it was also ~ 1 dex 
overabundant. 
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Fig. 7. As Fig. 6 for lines of Ni and s-process elements. 
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Fig. 8. Simulations of the blue region (4398— 4640 A) of the spectrum of FGSge between 1965 and 1972 assuming (a) solar abundances and 
(b) abundances from Langer et al. (1974). Each panel shows simulations with Tea decreasing from 8 000 K (Sp A8I: top) to 5 500 K (Sp GOI: 
bottom) in steps of 500 K, and has been plotted for comparison with Fig. 1 of Langer et al. (1974). The models have been degraded to an 
instrumental resolution of 0.25 A. Note the progressive enhancement with reduced T^s of lines due to s-process elements compared with lines 
of iron-peak elements. 



The Y and Zr abundances are a mystery. Without dredge- 
up there should be no change in these elements. With dredge-up 
there should be large changes in many other elements. Blocker 
& Schonberner (1997) assumed that some elements could have 
been locked in grains, but that orily moves the problem to one 
of how the grains would be dissociated as the star cools. 



5.6. Spectral evolution in the 1970's 

Cohen & Phillips (1980) give equivalent widths for a limited 
number of strong lines measured from echelle spectra obtained 
on 1975 June 24, 1976 July 7, and 1978 June 17. Using the 
new grid of model atmospheres, elemental abundances were 
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Fig. 9. Elemental abundance excess as a function of Tefi derived from the equivalent widths of Cohen & Phillips (1980). Solid and dashed lines 
refer to vt = 10 and 5kms~'^ respectively. Abundances are averaged over the spectra obtained in 1975, 1976 and 1978 and computed for a 
normal hydrogen abundance. The ordinate [Z] refers to log abundance (by number) normalized to the solar value. The + symbol represents 
solar abundance at the Teff appropriate for FG Sge in this epoch. 



computed as a function of T^a, Vt and tih- Figure 9 shows the 
results for normal nn only. In line with arguments presented 
elsewhere, normal uh, vt = lOkms"^ and Jeff = 5 000K 
are the most appropriate to use for these data. The main conclu- 
sions are summarised below. There are uncertainties associated 
with the equivalent width measurements and the atomic data 
used in the line formation calculations. These are illustrated by 
a spread of 2 dex in the results from five Y II lines. Currently the 
only Y I line and two Fe I lines measured by Cohen & Phillips 
(1980) give abundances which are discrepant from other hues 
due to the same atoms by a further 2 dex. There is little or no ev- 
idence for a change in composition during the interval spanned 
by these data. 

The individual line analyses are supported by examination 
of one echelle order of the 1975 spectrum of FGSge (Cohen 
& Phillips 1980, Fig. 3). Model spectra have been computed 
with the abundances of s-process elements and lanthanides Sr, 
Y, Zr, Ba, La, Ce, Pr, Nd, Sm and Eu enhanced by 0, 1 and 2 
dex in turn for a range of Tes, tih and Vt. Figure 10 shows 
the models for normal hydrogen, Tgs = 5 500 K, log g ~ 0.4 
and Vt = lOkms"^. A comparison with Fig. 3 of Cohen & 
Phillips (1980) demonstrates that all of the principal features 
in the observed spectrum can be reproduced with an excess of 
~ 1 dex in the s-process elements except, of course, for the 
interstellar component of the NaD line. 



The conclusions from the line analyses in the latter half of 
the 1970's include: 

- Na, Mg and Fe were approximately normal; 

- Y and Zr were overabundant by 1-2 dex, subject to a large 
line-to-line scatter in Y; 

- Ba, Ce, Pr and Nd were overabundant by ~ 1 dex; 

- La was overabundant by ~ 2 dex. 



5.7. Teff in the 1980's 

There are three published views of the run of T^s after 1980. 
Visual photometry provided either a ~ 4 500 K or ~ 5 500 K 
solution (van Genderen & Gautschy 1995), whilst the ultravi- 
olet fluxes observed with lUE led to a still higher ~ 6 500 K 
solution (Montesinos et al. 1990). The latter are based on a very 
small fraction of the total energy distribution of the star. This 
may be verified by combining the UV data with the visual and 
infrared data presented by the same authors (Table 2), and by 
comparing these with theoretical spectral energy distributions 
(SED's,Fig. 11). 

Best fit model energy distributions have been obtained us- 
ing a minimization procedure (FFIT, Jeffery 2001). We have 
used SED's from the models computed for this paper and have 
assumed Eb-v = 0.40 (van Genderen & Gautschy 1995) and 
a standard Galactic extinction law (Seaton 1979). Only lUE 
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Table 2. Effective temperatures and angular radii for FG Sge from spectral energy distributions. 



Date 


lUE 


V 


J 


H 


K 


L 


TeB 


xl 




image 


(FES) 










(kK) 


(/10~" rad) 



1982 June 


LWR 13399 


9.4 


7.0 


6.8 


6.6 


6.6 


5.45±0.08 


5.15±0.05 


0.47 


1986 October 


LWP 09356 


9.1 


7.0 


6.6 


6.4 


6.1 


5.42±0.05 


6.08±0.02 


3.05 


1987 July 


LWP 11269 


9.2 


7.0 


6.7 


6.5 


6.2 


5.39±0.05 


5.42±0.05 


0.91 


1987 Sept 


LWP 11585 


9.1 


7.0 


6.6 


6.5 


6.3 


5.30±0.04 


5.62±0.05 


0.57 


1987 October 


LWP 11820 


9.2 


7.1 


6.7 


6.5 


6.1 


5.28±0.04 


5.36±0.05 


0.84 
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Fig. 10. Synthetic spectra in the region of the NaD lines computed 
for comparison with Fig. 3 of Cohen & Phillips (1980). All spectra 
are calculated for a model atmosphere with Ted — 6 000 K, log g = 
0.5, Vt = lOkms^^ and convolved with a 0.4 A FWHM Gaussian. 
Three different compositions are assumed: (1) normal (dotted/red), (2) 
s-process elements enhanced by 1 dex (solid/black) and (3) s-process 
enhanced by 2 dex (dotted green). 
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Wavelength (A) / 10* 

Fig. 11. Model atmosphere fits to the observed spectral energy distri- 
bution of FG Sge in 1986 October. The main panel shows the UV flux 
distribution from lUE image LWP 09356 (histogram, lower left), the 
visual and infrared fluxes (Montesinos et al. 1990), and a series of 
theoretical models for Tctt — 5 000 (dotted, cyan), 5 500 (solid, blue) 
and 7 000K (dashed, green). The first two are normalised to the ob- 
servations at 5500 A, the latter is the optimum fit to the UV data. The 
bold (red) line is the best fit solution from FFIT cited in Table 2. The 
inset panel shows an enlargement of the same figure in the lUE range. 



data between 2400 and 3 100 A were used, larger data errors 
at A > 3100 A substantially degraded the fit quality. We could 
not differentiate qualitatively between models having 1%, 10%, 
and 90% hydrogen (by number); this introduces a systematic 
uncertainty in Tcs of a few tens of degrees, comparable with 
the formal fitting errors quoted in Table 2. Random errors in 
Tcs due to other factors (primarily data quality) are '--^ ±100 K, 
and systematic errors due to other properties of the models may 
also be ~ 100 K. 

We conclude that the high Tcs photometric solution of van 
Genderen & Gautschy (1995) is consistent with the ultraviolet 
and infrared photometry. We suspect that the measurements of 
Montesinos et al. (1990) were systematically in error due to 
the inclusion of noisy lUE data at A > 3100 A and since they 
are completely inconsistent with the contemporary VHJKL 
photometry. 



5.8. Did s-process abundances continue to 
increase into tlie 1990's? 

The most recent analysis of the FG Sge photosphere before it 
was obscured by circumstellar dust finds a continued increase 
in s-process and rare-earth elements (Gonzalez et al. 1998). 
Section 3.6 discusses how this conclusion was based on a mea- 
surement of Toff = 6 500 K which both contradicts the total 
flux distribution (Section 5.7) and implies a rapid contraction 
between 1992 and 1994. 

Prof. G. Gonzalez has kindly provided a copy of his 
1994 September 21 spectrum for comparison with our mod- 
els. While, we have not had the resources to carry out a com- 
plete fine analysis of this spectrum, we have been able to test 
(and revise) various hypotheses. As before, we have assumed 
logL/Af — 4.0 and have computed theoretical spectra for 
Toff in the range 4 500 to 7 000 K. We have examined mod- 
els with Vt = 5kms~^ and 10 kms~^, and three hydrogen 
abundances, rin = 0.90, 0.10 and 0.01. We have also consid- 
ered models with elemental abundances for Sc, Ga, Sr, Y, Zr, 
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Fig. 12. Model spectra in the region of Ha for several values of T^b (labelled) with log L/M = 4.0. The top row shows models with normal 
hydrogen abundance (nu = 0.90), and with s-process abundances increasing from normal (heavy black) through +1 dex (red) to +3 dex 
(green). The second row shows hydrogen-deficent models (nu = 0.01), and with s-process abundances increasing from normal (heavy black) 
through +1 dex (red) to +2 dex (green). The dotted line is the 1994 October 21 spectrum from Gonzalez et al. (1998). The model spectra have 
been convolved with a 0.5 A FWHM instrumental broadening function. Lines identified are from Kipper & Kipper (1993). 



Mo, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Er, Yb, Lu, Hf 
and Pb increased by 0, 1,2 and 3 dex. 

Figure 12 is plotted for direct comparison with Fig. 14 of 
Gonzalez et al. (1998). It shows a subset of the above models 
together with the spectrum of 1994 Sept 21 around Ha. The 
wavelength region shown includes lines of several s-process 
elements (La, Nd, Pr and Sm) as well as other species (Si, V, 
Cr, Ti). 

With normal hydrogen (riH = 0.90), T^s = 6 500K, 
log L/M — 4.0, and a 3 dex enhancement of s-process el- 
ements, the model does a reasonable job of reproducing the 
observed spectrum, except that the s-process lines are some- 
what weak and the region close to Ha, especially Sil, does 
not fit well. A better fit to the s-process lines is recovered with 
Teff = 5 500K,logL/Af = 4.0, and [s/Fe] = +3 redwai'd of 
Ha. However the "continuum" and the region to the blue are 
not reproduced. 

Models with Tcff = 5 500 K match the measured lUE 
fluxes. 

With reduced hydrogen ( nn 0.01), T^s = 5 500 K and 
log L/M = 4.0, several features are already quite well re- 
produced with normal (solar) abundances. Indeed, assuming a 
near-normal abundance for silicon, it appears in this case that 



the observed strength of Sil is a good indicator of hydrogen- 
deficiency. Exceptions to the fit quality are at the red and 
blue extremes. The remaining features can be reproduced with 
1 < [s/Fc] < 2. These lower [s/Fe] values result entirely from 
the reduced hydrogen abundance, which lowers the background 
opacity and immediately strengthens all of the absorption lines 
(for weak lines, line strength varies as the ratio of line opacity 
to local continuous opacity). 

Unfortunately, we are not in a position to model the C2 
bands and hence to establish whether the carbon abundance is 
significantly different from normal. The mere fact, however, 
that C2 is apparent indicates with certainty that C/0 > 1 and 
that the carbon abundance must be somewhat enhanced (see 
Section 3.5). 

In summary, by adopting a more appropriate T^b and grav- 
ity the FG Sge spectrum in 1992 shows: 

- reduced hydrogen abundance, 

- aC/Oratio > 1, 

- modest (^1 — 2 dex) overabundances of La, Nd, Pr and 
Sm, 

- by extension, detections of 3 dex overabundances in other 
elements seem premature. 
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Table 3. Blueward evolution of a 0.625 M0 LTP model, including 
principal abundances written as N^i (adapted from Blocker 2001). 



TcS 


logs 


Age (yr) 


H 


He 


C 





6700 


0.55 





0.918 


0.079 


0.003 


0.000 


7000 


0.63 


128 


0.877 


0.119 


0.003 


0.001 


7900 


0.80 


161 


0.793 


0.182 


0.021 


0.005 


9000 


1.00 


195 


0.657 


0.279 


0.052 


0.013 


10000 


1.20 


279 


0.577 


0.335 


0.071 


0.017 


11200 


1.37 


374 


0.461 


0.416 


0.099 


0.024 


11900 


1.47 


418 


0.321 


0.513 


0.134 


0.032 


12300 


1.49 


507 


0.181 


0.607 


0.170 


0.041 


12600 


1.49 


788 


0.162 


0.620 


0.176 


0.042 


13800 


1.70 


1073 


0.162 


0.620 


0.176 


0.042 


15800 


1.94 


1098 


0.162 


0.620 


0.176 


0.042 



A complete reanalysis of this spectrum using appropriate 
model atmospheres would be rewarding. 

5.9. Summary 

Taken as a whole, this re-evaluation provides the following 
view of the evolution of FGSge between 1960 and 1995: 

- The high T^b scale (van Genderen & Gautschy 1995) based 
on visual photometry and pulsation period proxies is con- 
sistent with lUE and infrared photometry, and excludes a 
contraction/heating phase from the 1980s's onwards. 

- FG Sge appears to have been hydrogen-rich at the start of 
this interval, but hydrogen-poor by the end. 

- There is little evidence for an unusual carbon abundance 
before the appearance of C2 bands in the 1990's, which 
demand C/O > 1. However, the available data are sparse 
(Herbig & Boyarchuk 1968) and the C abundance in 1960 
is poorly defined. 

- In general, the abundances of s-process elements and rare- 
earths are consistently about 1 dex above the solar value. 
With one exception, the evidence for wholescale changes 
of up to 3 dex must be viewed with scepticism. A modest 
increase (up to ~ 1.5 dex) between 1975 and 1994 may 
be acceptable, representing a confidence level on our esti- 
mates. 

- The only credible exception is a 1 dex enhancement in Y 
and Zr between 1965 and 1972. This is a mystery without 
dredge-up. 

6. Evolution 

The evolution of FGSge is defined photometrically (van 
Genderen & Gautschy 1995) since about 1890 and spectro- 
scopically since 1960. An additional upper limit of Tes — 
50 000 K before expansion is provided by the associated plan- 
etary nebula (He 1-5: Hawley & MiUer 1978), taking into 
account the weakness of the observed He 11 emission and 
the recombination time scale. The nebula is reported by the 
same authors to have normal abundances, although a previ- 
ous report of anomalously high N/O coupled with Jeff ^ 



55 000 ... 75 000 K should be acknowledged (Harrington & 
Marionni 1976). These constraints may be discussed in terms 
of theoretical models for the evolution of FG Sge. 

6. 1. Late thermal pulse models 

The upper limit on Tes defined where the helium-shell flash 
occurred during FG Sge's contraction towards a white-dwarf 
configuration and indicates that this event happened along the 
horizontal part of its post-ABG evolution. This type of flash 
is normally called a 'late' thermal pulse (LTP) as opposed to 
a 'very late' thermal pulse (VLTP) which occurs on the de- 
scending part of the post-AGB track (cf. Blocker 2001). From 
a theoretical point of view the characteristic difference between 
both types of thermal pulse is that only for the VLTP does 
the pulse-driven convective shell penetrate into hydrogen-rich 
layers, thereby causing complete mixing and burning of the 
hydrogen-rich envelope matter (Schonberner 1979; Iben et al. 
1983; Iben 1984; Herwig et al. 1999). According to flie latter 
authors, the final surface abundance pattern following a VLTP 
should reflect virtually that of the intershell region, viz. a mix- 
ture of helium, carbon and oxygen, without hydrogen (see also 
Herwig 2001a). 

The speed of the redward looping depends on the char- 
acteristic thermal/gravitational time scale of all the matter 
above the helium shell and is a strong function of the rem- 
nant mass. Based on 'flashing' post-AGB models of different 
masses Schonberner & Blocker (1997) estimated for FGSge a 
mass very close to 0.6 Mq. 

It is important to note here that in the VLTP mixing and 
burning of the hydrogen-rich envelope occurs immediately dur- 
ing the flash and is already completed before the star turns back 
to the red-giant region. On the contrary, a central star experi- 
encing an LTP is not expected to show any signs of internal 
mixing/burning processes at its surface while it is expanding 
towards a 'bom-again' red giant configuration. This appears 
to be confirmed by the normal surface composition found for 
FG Sge after it attracted the attention of observers, i.e. at a time 
when the flash is already over. 

In any case, the observed enrichment of the s-process ele- 
ments need not be the result of the late thermal pulse (cf. Iben 
1984). The s-process elements are likely to be the result of 
the past pulses when FG Sge was still an AGB star. Indeed, 
several post-AGB stars are already known to be enriched in C 
and N, with C/O > 1, weak iron-group elements and s-process 
elements by - 1 - 2 dex (Reddy et al. 1999, 2002, Hrivnak 
& Reddy 2003). With very similar surface compositions, such 
stars are Ukely counterparts to FG Sge should they suffer a LTP 
sometime in the future. 

Though the standard evolutionary models of an LTP de- 
velop strong envelope convection while approaching low sur- 
face temperatures, the bottom of the convection zone does not 
penetrate into nuclearly processed layers, i.e. the 3rd dredge-up 
does not occur until after the star reaches its minimum temper- 
ature. 

This is iUustrated in Fig. 13 and Table 3 by means of 
the 0.625 bom-again model of Blocker (2001). Dredge- 
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Fig. 13.Evolutionary track of a 0.625 M© 
post-AGB model suffering from a late ther- 
mal pulse (LTP). Symbols along the track 
give the time elapsed in units of lO^y, 
with zero set to the minimum effective 
temperature after the heUum shell flash. 
Dredge-up commences shortly thereafter. 
Age — 3 500yrs corresponds to the begin- 
ning of the post-AGB evolution. The figure 
is from Blocker et al. (2001). 



up starts after the reddest point is reached, leading to a slow 
dilution of the hydrogen-rich envelope by helium/carbon-rich 
intershell matter, with a mean rate of ~ 10~^ M0/yr. Note 
that there occurs no burning of hydrogen. During the dredge- 
up the star contracts slowly towards T^y = 10 000 K, at a 
rate of 10 K/yr. The dredge-up is completed within about 800 
years and leads to a hydrogen-deficient stellar surface: X = 
0.05, Y = 0.45, Xc = 0.38, Xq = 0.12. The model has 
obviously changed into a [WC]-type central star. A very sim- 
ilar result has been obtained by Herwig (2001) for a 0.6 Mq 
post-AGB model. 

Therefore the dilemma posed by FG Sge is that, according 
to the associated planetary nebula, it must have experienced 
a LTP. However, LTP models predict that surface hydrogen 
should not be depleted until a deep convection zone is estab- 
lished, and then only ~ 10* y after the pulse occurs. In the 
case of FG Sge, this hydrogen-depletion has occurred before 
the star reaches its minimum temperature, but substantially af- 
ter the time at which it is predicted to occur following a VLTP. 

6.2. Double-loop models 

An alternative and interesting model has been proposed by 
Lawlor & MacDonald (2003). In their calculations of the evo- 
lution of a post-AGB star following a 'very late thermal pulse', 
"models with low convective mixing efficiency first evolve 
quickly to the AGB, return to the blue, and then evolve more 
slowly back to the AGB for a second time before finally return- 
ing to the white dwarf cooling track". The attraction of this 
model is that it appears to provide a very satisfactory explana- 
tion for the rapid redward evolution of V4334 Sgr (Duerbeck et 



al. 2002), the historical appearance and subsequent contraction 

of V605 Aql (Clayton & de Marco 1997) and the comparatively 
leisurely expansion of FG Sge as different phases of the same 
evolutionary track. 

According to this model, FG Sge would be in its second 
return to the red, but the weaknesses of the model are several. It 
predicts that the surface of FGSge should be highly-processed 
(i.e. hydrogen-deficient and carbon-rich) throughout its recent 
history which, we have shown, it was not. It also states that the 
star would have been a white dwarf, with T^.^ > 100 000 K, 
in the comparatively recent past; too hot to comply with the 
nebula diagnostics (cf. Sect. 6), and possibly too old for the 
nebula to have survived (which it has). 

Very recently. Miller Bertolami et al. (2006) have made 
an independent calculation of evolution following a very late 
thermal pulse, and again find a double-loop structure with hy- 
drogen and then helium providing two phases of expansion. 
The interepretation of FG Sge within this framework suffers 
the same difficulties as for the Lawlor & MacDonald (2003) 
models. 

7. Conclusions 

We have made an extensive review of the chemical evolution of 
the surface of FG Sge between approximately 1960 and 1990, 
an epoch during which there have been published reports of 
dramatic changes in the abundances of many elements includ- 
ing particularly those created by s-process nucleosynthesis in 
the stellar interior. In several cases, these various reports have 
been shown to be seriously inconsistent with one another and 
with stellar evolution theory. Although it is more usual to recast 
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the theory when observation and theory come into conflict, the 
observational inconsistencies are sufficiently strong that they 
must be examined closely. 

Although we have not had direct access to the majority of 
the raw spectrograms, we have been able to assess the reliabil- 
ity of the majority of key studies, particularly those by Herbig 
& Boyarchuk (1968), Langer et al. (1974), Cohen & Phillips 
(1980), and Gonzalez et al. (1998). Although the first three of 
these report some changes in the abundances of carbon and s- 
process elements, we find that the data presented are almost 
entirely consistent with a constant surface composition typi- 
cal of other post-AGB stars, i.e. with C/O > 1 and s-process 
abundances enhanced by ~ 1-2 dex. The exception that the 
earliest data for Y and Zr indicate an increase from normal 
abundance to ^ + 1 dex remains a mystery. The very dramatic 
increases in s-process abundances reported in the most recent 
study (Gonzalez et al. 1998) are demonstrably wrong. They 
are due to the adoption of a value for T^s which is 1 000 K 
too high. A reassessment of the calibration of the photomet- 
ric evolution from ultraviolet data by Montesinos et al. (1990) 
demonstrates that the "high" Ten calibration by van Genderen 
& Gautschy (1995) is most consistent with all of the photomet- 
ric and spectroscopic data. Since FG Sge has now maintained 
its current Tcs of around 5 500 Kfor over a decade, this may 
represent the cool limit in its red ward evolution. 

Supporting Iben (1984), we have argued that the current 
surface of FG Sge does not show evidence of freshly dredged- 
up material in the form of carbon or s-processed elements. 
Contrary to Schonbemer & Jeffery (2003), we find that sur- 
face hydrogen has been depleted in the recent past (~ 40 y). In 
the 1960's it had an apparently normal H content (cf. Fig. 6), 
whereas in 1994 there is evidence of H-depletion via a reduced 
background opacity (cf. Fig. 12); a conclusion based on a par- 
tial analysis of the high-resolution spectrum of Gonzalez et al. 
(1998) not earlier at our disposal. We do not see evidence of 
the very prompt mixing suggested by recent VLTP theory. 

Stellar evolution theory predicts that, although a convective 
envelope does form as the star cools following a late thermal 
pulse, a deep convective envelope capable of dredging up pro- 
cessed material will not form until some time after the star has 
completed its redward evolution. Therefore, the observed hy- 
drogen depletion, suggesting prompt envelope mixing follow- 
ing a VLTP, contradicts the age derived from the associated 
planetary nebula, which suggests an LTP The question is then 
how well the physics of the LTP and VLTP are understood, and 
whether the demarcation between them is as concrete as current 
models suggest. 

FGSge is one of those very rare objects which demon- 
strates stellar evolution "while-you-wait". As a very evolved 
star, its surface already shows some material which has 
been created by nucleosynthesis in the deep interior. Astero- 
archaeology is the science of sifting through such material for 
clues to the past evolution of the star. In this case, most of the 
material was brought to the surface when the star had a deep 
convective envelope on the asymptotic giant branch. However, 
it seems that FG Sge has recently entered a new period of exca- 
vation, which has recently brought fresh helium to the surface. 
Other elements will be added to this within the next hundred 



or so years, when it will be possible to observe fresh material 
being brought to the surface as the star begins to dig up its own 
past. Continuing observations will provide key tests for stellar 
evolution models, particularly for models of late and very late 
thermal pulses. 
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Appendix A: Herbig & Boyarcliuli abundance 
measurements 

In order to make the abundance measurements for FGSge 
presented by Herbig & Boyarchuk (1968) more accessible, 
they are presented in Table A.l in the more conventional 
form of the logarithm of the abundance by number relative to 
a Cyg. Adopting a recent measurement of abundances in a Cyg 
(Albayrak 2000), the mean Herbig & Boyarchuk (1968) abun- 
dances in FG Sge are also shown relative to the Sun. Note that 
Albayrak (2000) gives the Cii abundance in a Cyg as 8.94 from 
one line only. The Ci abundance is 8.21. Since the Herbig & 
Boyarchuk measurements were based on Cll lines, we have 
used the Cll value for a Cyg in Table A.l. The mean over- 
abundance of C relative to the Sun between 1960 and 1965 in 
FG Sge would be H-0.4 dex rather than H-l.l dex if the Cl value 
were adopted for a Cyg. 
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Table A.l. Elemental abundances in FGSge: a) relative to aCyg (Herbig & Boyarchuk 1968) (cols 2 - 6), b) mean values (1960 - 1965, 
excluding 1960 Sept.) normalized such that the combined abundances of Si, Ti, Cr and Fe are the same in FG Sge and a Cyg (Albayrak 2000) 
(col 8), c) the same mean values relative to a Cyg (col 7) and the Sun (col 9), and d) elemental abundances for a Cyg (Albayrak 2000) and the 
Sun (cols 10 - 11). Abundances [X] are given as logarithm of relative abundance by number normalized to H = 12.00 
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